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We report on the observation of single top quark production by the DO collaboration using a 
dataset of 2.3 fb _1 collected at the Fermilab Tevatron pp collider. Several multivariate techniques 
are combined to separate the single top signal from backgrounds. The measured single top cross 
section is O t b+tqb = 3.94 ±0.88 pb. The probability to measure a cross section at this value or 
higher in the absence of signal is 2.9 x 10~ 7 , corresponding to a 5.0 standard deviation signifi- 
cance for the presence of signal. The lower limit at the 95% C.L. on the CKM matrix element V t b 
is \V,b\ > 0.78. A separate measurement of the f-channel cross section gives O tq b = 3.14+qjq pb. 
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Top quarks were first observed as pairs produced via the strong interaction at the Fermilab 
Tevatron Collider in 1995 ^ §]. Single top quark production proceeds via the weak interaction 
and its production cross section provides a direct measurement of the the quark mixing matrix 
element \V t h\ ®. It also serves as a probe of the Wtb coupling [Q, |[ |6[ [7|, ||] and is sensitive to 
several models of new physics [gp. 

In 2007, DO presented the first evidence for single top quark production and the first direct 
measurement of \V t b \ using 0.9 fb _1 of Tevatron data at a center-of-mass energy of 1.96 TeV [ |I0| , 
11]. Recently, the CDF collaboration has also presented such evidence in 2.2 fb~ ! of data [O]. 
Here we describe the observation of a single top quark signal in 2.3 fb _1 of data [13]. The CDF 
collaboration has also reported observation of single top quark production [jl4j]. 

Single top quark production proceeds via the s-channel production and decay of a virtual 
W boson (tb) and the ^-channel exchange of a virtual W boson {tqb). We search for both of these 
processes at once as well as the ^-channel process alone. The sum of their predicted cross sections 
is 3.46 ± 0. 18 pb [0] for a top quark mass m, = 170 GeV. 

We select events collected with the DO detector [jl^] containing one isolated lepton (electron 
or muon), missing transverse energy, and two, three, or four jets, at least one of which is ^-tagged. 
We separate the analysis into 24 channels by lepton type, jet and ft-tag multiplicity, and run period. 
We model the signal using the CompHEP-based next-to-leading order (NLO) Monte Carlo (MC) 
event generator SINGLETOP JT7|]. The main backgrounds to the single top final state signature are 



W+jets and Z+jets production, as well as tt, all of which are modeled using ALPGEN [18p. A smaller 
background is from multijet events which are modeled using data samples. Other backgrounds are 
from diboson production, modeled using PYTHIA [|19|]. 

Systematic uncertainties arise mainly from the jet energy scale corrections and ft-tag model- 
ing, with smaller contributions from MC statistics, correction for jet-flavor composition in IV+jets 
events, and from the IV+jets, multijets, and tt normalizations. The total uncertainty on the back- 
ground is (8-16)% depending on the analysis channel, and we take both normalization and shape 
effects into account. 

We select 4519 events with a background expectation of 4428 events and a single top signal 
expectation of 223 events. Since the signal is small compared to the overwhelming background, we 
apply three separate multivariate analysis techniques based on boosted decision trees (BDT) [20, 
2ll [22|], Bayesian neural networks (BNN) pf, M\, and the matrix element (ME) method |M Eflft 



to extract the single top signal. We then combine these in a combination BNN. 
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Figure 1: Combination discriminant in two cross-check samples. 
The BDT analysis uses a common set of 64 discriminating variables for all analysis channels, 



2 



DO Single Top Observation 



Reinhard Schwienhorst 



but optimizes the filters separately in each channel. The BNN analysis selects 18-28 discriminating 
variables in each channel. The ME analysis uses only 2-jet and 3-jet events and splits the analysis 
into I0W-//7- and high-Z/y regions at Hj = 175 GeV. All three analyses transform their output distri- 
butions to ensure that every bin contains sufficient background statistics. We verify the agreement 
between data and background model for each multivariate method separate cross-check samples: 
a W+jets dominated sample and a ff dominated sample. The combination discriminant output for 
these two samples is shown in Fig. [[]. Fig. i5] shows the combination discriminant output together 
with one of the discriminating variables for events in the signal region. The measured cross section 
is Gtb+tqb = 3.94 ±0.88 pb. The measurement has a p-value. of 2.6 x 10~ 7 , corresponding to a sig- 
nificance of 5.0a. We use the cross section measurement to determine the Bayesian posterior for 
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Figure 2: Combination discriminant for all analysis channels combined (a), only the signal region (b), and 
the pseudorapidity of the light quark jet multiplied by the lepton charge for events in the signal region (c). 
The signal is normalized to the measured cross section. 

\V t b\ 2 in the interval [0,1] and extract a limit of \V t b\ > 0.78 at 95% confidence level. 
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Figure 3: Ranked f-channel discriminant (a), only the signal region (b), and posterior probability density 
for f-channel and .s-channel single top quark production (c). Also shown are the measured cross section, SM 
expectation, and several representative new physics scenarios [0, 53]. 

In a separate analysis we train the same multivariate methods using only f-channel single top 



events as signal and including SM s-channel in the background [27]. The f-channel discriminant 
output is shown in Figs. |3|(a) and ||(b). We then determine the posterior as a function of both s- 
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channel and f-channel as shown in Fig. |](c). From this posterior we obtain the ^-channel cross 
section by integrating over the s-channel axis and find o tq b = 3.14^Q gQ pb. We similarly extract the 
s-channel cross section as o t b = 1 -05 db 0. 8 1 pb by integrating over the £-channel axis. The observed 
p-value for the £-channel measurement is 8.0 x 10~ 7 , corresponding to a Gaussian significance of 
4.8a. 
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